Onset of Couette flow in a confined colloidal glass 
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Abstract 

A colloidal glass confined between rough walls is investigated using non-equilibrium molecular 
dynamics computer simulations. The onset of flow, under the imposition of a uniform shear stress, 
is studied. When the imposed stress is gradually decreased, the timescale for the onset of steady 
flow diverges. Near this yield-stress regime, persistent creep is observed. This regime of slow 
deformation is shown to be associated with long-lived dynamical heterogeneities, which are strongly 
influenced by the initial state. 
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A defining mechanical property of amorplious solids is the existence of a yield stress; these 
materials yield only when the applied stress exceeds this threshold In experiments, the 
material's response near yielding is commonly probed by using stress as a control parameter 
[2!. Such measurements report about the existence of a creeping regime, where the material 
deforms very slowly. Quite often, in such situations, steady flow is not observed in^soft 
amorphous materials (like colloids, gels, emulsions etc.) within experimental durations 
This is related to the recent recognition of diverging timescales for the onset of flow 
0|. Beyond these macroscopic measurements, little is known about the spatio-temporal 
characteristics of the local dynamics around yielding, which would allow for an improved 
understanding of the micro-mechanisms at play. 

An emerging scenario for the micro-dynamics is that when the material is sheared, struc- 
tural rearrangements occur locally js iO] which subsequently trigger more events in the 
neighbourhood and this correlated process initiates and sustains the flow. Kinetic elasto- 
plastic models 11| based on such a picture of fluidization have been able to explain various 
steady-state flow properties in conflned systems 12h14|. Further, it has been conjectured 
that the triggering process results in an avalanche-like behaviour, as evidenced via simu- 



lations 
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16|. Here, we focus on how such correlations are expressed when a quiescent 
glass is subjected to an external stress. The recent measurements near creep jsl seem to 
indicate that flow occurs via the formation of transient shearbands; however, whether these 

are consequences of 



flow inhomogeneities, which are also observed in steady flow |17l-20 



the above-mentioned correlated processes remains to be resolved [18|, |21|, |22|. Also, recent 
stress-controlled granular experiments 23|, while observing such shearbands, report the sig- 
niflcance of local plastic events leading to avalanches and eventual flow 2J] . So the question 
is whether these features are generic to the fluidization of amorphous systems under an 
imposed stress. 

In this Letter, we report a simulational study of the onset of flow in a conflned colloidal 
glass, under an imposed uniform shear stress. We use a geometry that mimics the typical 
planar Couette setup in stress-controlled experiments of soft materials. We show that indeed 
the timescale for the onset of steady flow diverges with decreasing stress, along with the 



emergence of nonlinear creep which has a power law form (similar to a range of materials 
25|-|28|). Moreover, by following the local dynamics, we observe that the mobility is 



m 



spatio-temporally heterogeneous and the mobile regions take the shape of shear-band-like 
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structures. Thus, our work provides evidence that in the presence of an external shear 
stress, the fluidization of the glass near yielding is extremely slow, along with the presence 
of long-lived spatial heterogeneities. 

For our simulations, the model colloidal system that we consider is a 50 : 50 binary 
Yukawa fluid 
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30| . Thus, interactions between the particles are given by the pair potential 



= Eapdajsexpl—K^r — dai3)]/r, with r the distance between a particle of type a and 
one of type (3 {a,P = A,B). For model parameters, see 3l|- Our MD simulations have 
been done for samples having the dimensions = 26.66(is, Ly = 53.32ds, = 13.33ds, 
consisting of = 12800 particles. We work in the NVT ensemble and the temperature 
is controlled by using a Lowe thermostat 32|. We equilibrate the system, using periodic 
boundary conditions, at a high temperature of T = 0.2 and then m = 24 independent 
configurations sampled at T = 0.2 were instantaneously quenched to T = 0.05 (which is 
below the mode-coupling critical temperature of = 0.14 |29|). Subsequently, at T = 0.05, 
each of these m configurations are aged for durations of t^gc = 10^, 10^, 10^. When each 
configuration reaches a certain t^gc, we freeze the particles at < y < 2ds and Ly — 2ds < y < 
Ly to prepare glassy states confined between rough walls. Then, these confined samples are 
sheared by pulling the top plate in +x direction with a fixed force Fq (similar to rheometers) 
which imposes a constant shear stress of o"o = Fq/L^Lz] we study the response of the 



confined glass to such an imposed field. Alternatively, a constant shear-rate (70) can be 
imposed by pulling the top plate at a constant velocity vq = 'jo{Ly — Adg) |20|]; we used such 
a setup to estimate the dynamical yield stress (ad) of the confined glass at T = 0.05 |34| . 
obtaining a value of ad = 0.0857. 

Similar to experimental creep measurements, we study the response of the glass by mon- 
itoring the strain experienced at the top wall 7w(^)- This is done by recording, for each of 
the m samples, the velocity of the wall v^{t) as a function of time, and then the average 
effective strain at the wall is obtained: 7w(^) = (7a(^))e, where 'jait) = jQ[vu){t) / w]dt is the 
wall strain computed for a-th sample and ()e is an average over the ensemble of m trajecto- 
ries. In Fig. [11(a), we show how 7w(^) evolves for a wide range of imposed stress ctq- These 
curves show several regimes - (i) in all cases, at early times, the strain increases initially 
and then there is an oscillatory part, corresponding to the regime when the stress builds up 
inside the confined sample - the oscillations occur due to the interplay of the imposed stress 
at the wall and the restoring force of the deformed glass; (ii) after this, for large values of 
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FIG. 1: (a) At an age of tage = 10^, evolution of strain at wall, 7w(i) for a range of imposed stress 
(Tq. (b) Variation of timescale for onset of steady flow, tq, as a function of ctq. The dashed black line 
corresponds to a fit to the data with A/{ao — o"s)^, with as = 0.0848 and /3 = 2.285. (c) Variation 
of strain fluctations across trajectories, 6j^{t) for different do. The dashed line corresponds to a 
power law fit 6j„{t) ~ t^/^^ 



(To, one quickly sees an asymptotic linear regime in 7w(^), which corresponds to a steady 
flow at a fixed shear-rate; (iii) as ctq is decreased, an intermediate regime in 7w(i) emerges, 
which for small enough (Tq has a power-law behaviour - for the creep at ctq = 0-80 (which is 
below the dynamic yield stress 0"^), we obtain a power-law exponent of 0.25. For imposed 
stresses larger than ad, one eventually observes the asymptotic linear regime. Now, one can 
define a time-scale for onset of flow, tq, as the time required for the system to reach a strain 
of 7w(to) = 1; using this definition, we calculate tq for each ctq and plot it in Fig. [D^b). One 
can see that the timescales increase with decreasing stress and the data can be fitted with 
the function A/{ao — o"s)^, with = 0.0848 (~ ad, the estimated dynamical yield stress) 
and P = 2.285 (which is similar to experimental observations 5|, l7|). 

In supercooled liquids, diverging timescales are seen to be associated with increasing 
dynamical heterogeneities. For the onset of flow in glasses, we explore this possibility by 
studying the fluctuations in response within the ensemble of m samples having the same 
age; this is quantified by calculating 6'y^{t) = {la{t))e/ {laif))1 — 1- In Fig- II^c), we plot 
57w(t) for different ctq. We see that for large values of do, the function 57w(t) increases with 
time, has a maximum and then decays at long times. The location of the maximum shifts 
to larger timescales with decreasing ctq, similar to the trend for tq. As we approach the 
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yielding regime, the location of this maximum goes beyond our time window of observation. 
Interestingly, the initial increase in 57„(t) has a power-law behaviour, which is identical for 
all (Jo, as demonstrated in Fig. HJ^c). This is similar to what was earlier observed in creep 



flow in experiments using paper samples 
power law in the fluctuations. 




- a power-law in average strain, as well as a 
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FIG. 2: (a) Wall-strain 7w(t) for samples having different ages (t^g,, = 10^10^10^) for different 
imposed stresses cjo = 0.115,0.100,0.092. (b) Corresponding effective strain-rates, 'y^{t)/t, at the 
wall, (c) 7w,(t), for (Jq = 0.080 for the three ages, revealing identical power-law creep (with exponent 
of 0.25, indicated with dashed lines) in this sub-yielding regime. 

Aging is a characteristic property of glass-forming systems around the glass-transition 
temperature Tg. Below this temperature, the transient response to a deformation does 
depend on the history of preparation 33|,l35[. This has also been observed in the recent creep 
experiment on colloidal glass [3] , where the deformation of samples having different ages was 
studied. In our simulations, the ages that we can explore are relatively smaller compared 
to experiments; despite that, we explore whether similar aging-dependence is also observed 
in our case. In Fig. |2]^a), for different values of imposed stress (Tq = 0.092,0.100,0.115, we 
plot the wall-strain 7w(t) for samples having different ages, viz. tagc = 10'^, 10^, 10^. The 
corresponding effective strain-rates at the wall, 7w(t)/t, are shown in Fig. [2|^b). As observed 
in the experiments, for (Tq = 0.100,0.115, we recover the same asymptotic regime for 7w(^), 
i.e. the final shear-rate is the same for all the aged samples, which is expected. However, 
the transient regime depends upon the age - the larger the age, the later it reaches the 
asymptotic linear regime with variation in the shape of '~fw{t)/t when the system "breaks" 
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FIG. 3: For a = 0.080, spatial maps of average transverse displacements along a single trajectory 
(tagc = 10''). Measurements done at t = 10790,25730,92960,251905 {from left to right), with the 
colorbar giving the scale of local mobilities. 

into flow. In fact, the curvature in 7„(t)/t can also be different, which we observe for 
(To = 0.10 (see Fig. Mjo)) - while the curve for tage = 10^ has a pronounced minimum, it is 
a monotonically decreasing function for tage = 10'^. This implies that the intermediate flow 
states could be different for samples having different ages. For the case of ctq = 0.092, 7w(t) 
for tage = 10^ does not reach the steady-state within our window of observation, unlike the 
cases of tagc = 10^, 10^. In Fig. EJ^c), we study the situation for the case when o"o < (Td, for 
all the three different ages, at ctq = 0.080, 7w(^) has a power-law creep, with an exponent 
of 0.25. This indicates that the origin of the power-law behaviour is due to processes which 
are age- independent. 

Until now, we have measured the macroscopic response by monitoring the motion of 
the top wall. We observed, with decreasing stress, the appearance of creeping dynamics. 
We now want to see how this slow deformation shapes up at a more local scale within the 
sample; whether the local response is uniform or heterogeneous and whether any spatial 
structures are formed as the system yields. In order to study this, we consider non-affine 
motions occurring due to local structural rearrangements. To this end, maps of transverse 
displacements are computed in the following manner. Before we apply the external stress 
{t = 0), we divide the xy plane of the simulation box into small square cells (of length 
Ic = l-lds) and identify the particles in each cell. Next, after time t, we determine the 
transverse displacement of each particle Ayi{t) = \yi{t) — yi{0)\ and then, to construct the 
spatial maps, we calculate for each cell the local mobility, Himif) = {^yi{i))iTn (with the 
average over all the particles in the cell {/m} at t = 0). 

In Fig. [3l for a single trajectory of a sample aged to tage = 10^, we show the displacement 
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FIG. 4: For samples at tage = 10^, {left to right), the first eight panels show displacement maps for 
different trajectories at a fixed time t = 42745 for do = 0.092, indicating the variation in response 
of different samples at same age. {Last panel) 6j„{t) for the corresponding ensemble, with the 
arrow indicating the time of measurement. 

maps as the system evolves in time, at an imposed stress of (Tq = 0.080 (i.e. where we 
observe the power-law creep). Note that the local non-affine dynamics is extremely slow; 
even at t = 10790 (first panel on the left), only a few faint spots of significant displacements 
{lJ^im{t) > 0.5ds) are seen. As time progresses, these spots expand into larger patches {t = 
25730). By t = 92960, we see that these patches form a band-like structure extending across 
the length of the box and this structure remains persistent till the end of our observation {t = 
251905). In the rest of the system, even at such long timescales, the mobility is comparably 
negligible. Thus, during creep fiow, while initial mobile spots seem to initiate activities in the 
neighbouring region along the fiow direction, similar mobilities in the transverse direction 
seem to be impeded, which makes these localized structures persistent. For larger stresses 
where eventually steady-state fiow is obtained, such spatial heterogeneities also occur; but 
they are temporally short-lived, with the mobile regions quickly building up transversely 
across the sample until the entire system is fiuidized. Hence, the divergence of timescales for 
the onset of steady-state, with decreasing (Tq, is probably associated with a slow emergence 
of regions of mobility as well as the impediment of propagation of fiow in the transverse 
direction. 

The displacement maps can now be used to clarify the origin of the temporal behaviour 
of the fiuctuations in wall strain 6'y^{t). For samples aged to tagc = 10^, we study the 
dynamics for the independent trajectories evolving for o"o = 0.092. In Fig. Hl^a), eight maps 
from this ensemble, calculated at t = 42745, are shown; the time of measurement is marked 
by an arrow in Fig. Hl^b), which shows 57w(^) for the entire ensemble. From map-to-map, 
one clearly sees that the local mobility is spatially different, with the occurrence of varying 
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degrees of localization (or even the absence of it). Thus, this implies that even though all 
these confined samples have the same age, the local response to an imposed stress differs from 
one to the other, depending upon the initial state. The varying degree of fluidization, within 
an ensemble, indicates a distribution of timescales for complete fluidization of each initial 
state; and this variation is captured within 6^^{t). Therefore, the increase in S'~f^{t) with 
decreasing o"o implies that this distribution of timescales becomes broader as we approach 
the yielding regime. 



-'— h 

2 

3 

_ 4 


1 ' 1 


1 ' f 

(a) 


5 

- 6 

7 

8 

9 














1 , 1 


1 , r 


10 


20 30 


40 50 




FIG. 5: (a) Velocity profiles during the onset of flow at ctq = 0.10 for a sample with age tage = 10^, 
averaged for At = 498, starting from different time-origins Iq (labelled in sequence), (b) For the 
same trajectory, 7w(t) with the black squares marking the different to- 



The displacement maps, which track motion transverse to the flow direction, indicate the 
presence of heterogeneous dynamics; one can expect that the corresponding velocity profiles 
(which monitor the spatial dependence of the motion along the force direction) would also 
capture such heterogeneities. Recently, there has been focus on the possible occurrence of 
transient shear-bands [s], 3^ in velocity profiles during the onset of flow in soft glasses. For 
a single trajectory at tagc = 10^, we show such velocity profiles in Fig. [5]^a), measured during 
flow at ctq = 0-10 case where one attains a steady-state at long times). Since instantaneous 
profiles are too noisy, we average velocity profiles over short intervals in time (At = 498) 
starting at different points in time (to) after the imposition of the stress. In Fig. Mjo), we 
plot the evolution of the wall-strain 7w(^) for this trajectory and the different to are marked 
on the curve; they correspond to the timescales when the material "breaks" into flow. We 
see that, during these times, the velocity proflles fluctuate quite a lot, with signatures of 
intermittent flow (i.e. existence of little or no flow at some time instances, interspersed with 
burst of flow) as well as the existence of band-like profiles (co-existence of regions of no-fiow 
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with more mobile regions). These transient heterogeneities in the velocity profiles, while in 
agreement with the observation in gels jsl, seem to reflect a more complex motion in flow 
direction. 

Creep flow observed in crystalline materials, are often explained via the motion of dis- 
.ocat,™. H; Koweve. fo. the d-soMered .natena., such processes are absent. For soft 



materials, different rheological models like SGR 38|, phenomenological fluidity models 
and schematic models [3] do predict creep flows of various forms; but neither do they take 
into account any spatial information for building the model, nor do they provide predictions 



for the local dynamics. In that aspect, the non-local fluidity model seems promising; 
however, the temporal evolution of spatial patterns of local fluidity, as the material starts 
flowing, is yet to be worked out. Only recently, a more spatial version of the SGR model 



was used to show that creep flow is associated with shearbanded velocity profiles |40|. In 
this context, our simulations provide further insight into the possible ingredients for con- 
structing a more complete theoretical model which would correctly describe the dynamics 
near yielding. 

In conclusion, using numerical simulations, we have demonstrated that during the onset of 
Couette flow in glassy systems, creep (characterized by a power-law dependence of strain on 
time) is associated with strong spatial heterogeneities in the dynamics which take the form 
of long-lived shearbands. Moreover, for the confined system, the local response also depends 
on the initial structure, showing up as fluctuations within an ensemble. Thus, even for a 
simple planar Couette flow geometry, one observes a rich spatio-temporal response of the 
soft glass as it yields. Further studies are necessary to clarify how this transient behaviour 



connects to creep observed in uniaxial tensile studies of metallic [4l| and polymeric |42|-|44 1 
glasses. For soft amorphous materials, while more local measurements are necessary from 
experiments, simulations can also further explore yielding in different complex geometries. 
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